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Abstract. Many disordered crystalline materials show 
chemical short range order and relaxation of neighboring 
atoms. Local structural information can be obtained by ana- 
lyzing the atomic pair distribution function (PDF). The via- 
bility of reverse Monte Carlo (RMC) simulations to extract 
quantitative information about chemical ordering as well as 
displacements is investigated. The method has been applied 
to simulated PDFs of disordered structures showing chemi- 
cal disorder alone as well as in combination with displace- 
ments. 



Introduction 

The determination of crystal structures is an important part 
of chemistry, physics and of course crystallography. Con- 
ventional structure determination is based on the analysis of 
the intensities and positions of Bragg reflections which only 
allows the determination of the long range average structure 
of the crystal. For powder diffraction data this is now rou- 
tinely done using the Rietveld (Rietveld, 1969) method. It 



should be kept in mind that the analysis of Bragg scattering 
assumes a perfect long range periodicity of the crystal. How- 
ever, many materials are quite disordered and even more im- 
portant the key to a deeper understanding of their properties 
is the study of deviations from the ave rage structure o r the 
study of the local atomic arrangements ( Billinge, 1996 ). De- 
viations from the average structure result in the occurance 
of diffuse scattering which contains information about two- 
body interactions. It is interesting to remember that the in- 
formation that can be extracted from diffuse scattering is in- 
deed limited to two-bod y correlations. In an review article 
about diffuse scattering, Welberry and Butler (1994) show 
two chemically disordered structures giving rise to the same 
diffuse scattering pattern. Two similar structures are shown 
in Figure 0and one can easily see the different arrangement 
of the two atoms types shown as filled and empty circles. 
However, for both structure all two body correlations are 
zero, the only difference is a non-zero three -body term for 
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Fig. 1. (a) Structure with random defect distribution, (b) structure 
showing three-body correlations, but all two-body correlations are zero 
and (c) the resulting PDF of both structures given. The difference is 
magnified by a factor of 10 and plotted below the data sets. 



one of the structures, giving rise to the triangular arrange- 
ment seen in Figure |IJ In principle those two structures can 
not be distinguished using diffuse scattering or PDF mea- 
surements. The corresponding PDFs for the two structures 
are shown in Figure |l|c. However, in many cases three-body 
correlations might impose constrains on local displacements 
which then change the diffraction patter n or PDF. A deta iled 
discussion of this topic can be found in Welberry (1986 ). 

There are various ways to extract information about 
the defect structure from single crystal diffuse scattering. 
An overview of traditional approaches and the analysis of 
diff use scattering via computer simulations can be found 
in ( Welberry and Butler, 1994 ), further general informa- 
tion about disorder diffuse scattering can be found in nu- 
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merous review articles (|Welberry, 1985 


Tagodzinski, 1987; 


lagodzinski andFrey, 1993; prey, 1995, 


1997|; Welberry and 



Butler, 1995). However, the data collection, interpretation 
and analysis of diffuse scattering is a time consuming and 
in general difficult task. Information about the local envi- 
ronment of atoms in a given structure can be extracted us- 
ing extended X-ray fine absorption structure (XAFS) anal- 
ysis and related t echniques. A re cent review of those meth- 



ods was given by Yacoby (1999). The method to reveal the 



local structure of crystals used here is the analysis of the 
PDF. This approach is long known in the field of study- 
ing short range order in liquids and glasses but has recently 
been applied to crystalline materials dEgami, 1998]; T oby and 



Egami, 1992; [Billinge and Egami, 1993D - The PDF is ob 



tained from the powder diffraction data via a sine Fourier 
transform of the normalized total scattering intensity S(Q): 

= 4irr[p(r) - p ] 
2 f°° 

Q[S(Q)-l]sin{Qr)dQ, (1) 



G e (r) 



where p(r) is the microscopic pair density, po is the average 
number density and Q is the magnitude of the scattering vec- 
tor. For elastic scattering Q — An $m{6) / A with 28 being the 
scattering angle and A the wavelength of the radiation used. 
Since the total scattering S(Q) contains Bragg and diffuse 
scattering, the information about local arrangements is pre- 
served. Details about the determination of an experim ental 
PDF can be found e.g. in [Egami (1998| ); |Warren (1990| ) and 
are not discussed here. 

In many systems one is faced with chemical short range 
order and subsequent displacements of neighboring atoms. 
In previous work, the capability of the RMC technique to 
refine single crystal diffuse scattering of systems showing 
similar disorder was investigated (Proffen and Welberry, 
1997ag). However, obtaining a PDF from a powder diffrac- 
tion experiment is often simpler than carrying out a single 
crystal diffuse scattering measurement and allows one to 
e.g. collect data for many temperature points around a phase 
transition or explore a family of compounds more systemat- 
ically. The aim of this paper is to investigate the capability 
and possible limitations of reverse Monte Carlo (RMC) re- 
finement of PDF data to obtain information about chemical 
ordering and subsequent displacements. 



maximum distance), however each contribution has a weight 
corresponding to the scattering power of the two atoms in- 
volved. The PDF can be calculated from a structural model 
using the relation 



Gc(r) 



w s{r 



4irrp , (2) 



where the sum goes over all pairs of atoms i and j within 
the model crystal separated by ry . The scattering power of 
atom i is bi and (b) is the average scattering power of the 
sample. In the case of neutron scattering b$ is simply the 
scattering length, in the case of x-rays it is the atomic form 
factor evaluated at a user defined value of Q. The default 
value is Q = in which case bi is simply the number of 
electrons of atom i. 

Generally there are two different ways to account for dis- 
placements (either thermal or static) from the average po- 
sition. First one can use a large enough model containing 
the desired displacements and perform an ensemble aver- 
age. This can be achieved e.g. by displacing atoms in the 
model according to a given Debye-Waller factor. Alterna- 
tively one can convolute each delta-function in (||) with a 
Gaussian accounting for the displacements. The width try 
of the Gaussian is given by the anisotropic thermal factors 
Uim — (uiu m ) of atoms i and j. 

The study of a measured PDF ranges from a simple peak 
width analysis revealing information about correlated mo- 
tion (Jeong et al., 1999) to the full profile refinement of the 
PDF based on a structu ral model either using the Reverse 
Monte Carlo technique (Toby and Egami, 1992; Nield et al., 
1994) or least square regression ( |3illinge, 1998 ). The later 
can be carried out using the program PDFFIT (Proffen and 
Billinge, 1999). The reader might refer to this paper for de- 
tails on the calculation of PDFs and experimental factors like 
terminating the data at a given value of Qmax determined by 
the experiments. We will exclude the discussion of experi- 
mental influences on the PDF from this paper, since its scope 
is a principle investigation of the capability of RMC refine- 
ments of the PDF to obtain information about chemical or- 
dering and possibly accompanying displacements. All simu- 
lations and refinements presente d in this paper were carri ed 
out using the program DISCUS ( proffen and Neder, 1997 ). 



Simulation techniques 

In this section some details about the simulation techniques 
in general are given. In the first part the calculation of a PDF 
from a structural model is discussed, the following two parts 
describe the Monte Carlo (MC) and reverse Monte Carlo 
(RMC) simulation technique. Although this paper concen- 
trates on the RMC refinements of various simulated struc- 
tures showing occupational and displacive disorder, these 
defect structures themselves were created using the MC 
method. 

Calculating the PDF 

The PDF can be understood as a bond-length distribution 
between all pairs of atoms i and j within the crystal (up to a 



Monte Carlo (MC) 

The total energy of a model crystal is expressed as a function 
of random variables such as site occupancies or displace- 
ments from the average structure. The algorithm is based on 
the original Metropolis Monte Carlo technique (Metropo- 
lis et al., 1953). A site within the crystal is chosen at ran- 
dom and the associated variables are altered by some ran- 
dom amount. The energy difference AE of the configuration 
before and after the change is computed. The new configu- 
ration is accepted if the transition probability P given by 
equation || is greater than a random number h, chosen uni- 
formly in the range [0,1]. 

= cxp(-AE/kT) 

1 + exp(~AE/kT) K ' 
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In |] T is the temperature and k Boltzmann's constant. It 
should be noted that the value of the temperature T con- 
trols the proportion of accepted modifications which lead to 
a larger total energy. The process is repeated until the sys- 
tem reaches its equilibrium. In this paper we refer to a single 
MC (or RMC) step as 'move', whereas the number of moves 
necessary to visit every crystal site once on average will be 
called one 'cycle'. 

Next we briefly describe the actual energy expressions 
used to introduce occupational and displacement disorder 
using the MC algorithm. All simulations start from a given 
undistorted average crystal structure. Chemical ordering is 
described using binary random variables cr^ , with ov, = + 1 
to represent a Cu atom and <7y- = — 1 for a Au atom. The 
energy of interaction between sites used in the MC vacancy 
ordering scheme is of the form of 



(a) swap atoms 



En 



^2 <*ij [ H + ^IO-(IOO) + ^20-(ho> + ■■] (4) 



where O7 10 n) is the sum of all four nearest-neighbor vari- 
ables, (J(iio) the sum of all four next-nearest-neighbor vari- 
ables and so on. The sum in equation [| is over all crys- 
tal sites i and j. The interaction parameters H, J\ and J 2 
are initially unknown and a feedback mechanism is used to 
achieve the desir ed ordering of the vacancies(see W elberry 
and Butler, 1994). The displacements were modeled in a 
similar way using a Hamiltonian, where the atoms move in 
harmonic potentials (Hooke's law). 



E. 



disp 



id 



k(di 



A) 2 



(5) 



The sum is over all atoms i and all nearest-neighbor 
atoms j. Only nearest-neighbor terms are used in this mod- 
eling. The atom-atom distance is given by dij, the average 
distance is do, is the displacement factor and k a force 
constant. 



Reverse Monte Carlo (RMC) 

The RMC method was originally developed by McGreevy 
and Pusztai (1988). The principle is the same MC type al- 
gorithm described in the last section. However, rather than 
minimizing the total energy, the difference between ob- 
served and calculated data (in our case PDFs) is minimized. 
The algorithm works as follows: First, the PDF is calcu- 
lated from the chosen crystal starting configuration and a 
goodness-of-fit parameter \ 2 is computed. 



X 



j2 

i=l 



G c {n)f 



(6) 



The sum is over all measured data points r,, G e stands 
for the experimental and G c for the calculated PDF. The 
RMC simulation proceeds with the selection of a random 
site within the crystal. The system variables associated with 
this site, such as occupancy or displacement, are changed 
by a random amount, and then the model PDF and the 
goodness-of-fit parameter \ 2 316 recalculated. The change 
A\ 2 of the goodness-of-fit x 2 before and after the gener- 
ated move is computed. Every move which improves the fit 






• 


\ 








• 


• 


(b) shift atom 


O 


O 


O 


O 




>> 

o o 



(c) swap displacements 



9 





\ 











6 



6 


H0 





9 



Fig. 2. Illustration of the MC/RMC modes: (a) swap atoms, (b) shift 
atom and (c) swap displacements from average position. 



(A% 2 < 0) is accepted. 'Bad' moves worsening the agree- 
ment between the observed and calculated PDF are accepted 
with a probability of P = cxp(— A\ 2 /%)■ As the value of 
A% 2 is proportional to l/cr 2 , the value of a has an influence 
on the amount of 'bad' moves which will be accepted. Ob- 
viously there are two extremes: For very large values of a, 
the experimental data are ignored (x 2 ~ 0) and with very 
small values of a the fit ends up in the local minimum clos- 
est to the starting point, because there is a negligible prob- 
ability for 'bad' moves. The parameter a acts like the tem- 
perature T in 'normal' MC simulations. The RMC process 
is repeated until x 2 converges to its minimum. The result 
of a successful RMC refinement is one real space structure 
which is consistent with the observed PDF. In order to ex- 
clude chemically implausible resulting structures additional 
constrains, e.g. minimal allowed distances between atoms, 
may be introduced. 

The program used f or the RMC refinements p resented in 
this paper is DISCUS ( [Proffen and Neder, 1997] ). The pro- 
gram is capable of modeling occupational as well as dis- 
placement disorder and so far we have called each crystal 
modification simply RMC move. In practice we use a mode 
('switch-atoms') of simulation in which occupational dis- 
order is modeled by swapping two different randomly se- 
lected atoms (Fig. ga). This procedure forces the relative 
abundances of the different atoms within the crystal to be 
constant. It should be noted, that vacancies are treated as an 
additional atom type within the program DISCUS. The intro- 
duction of displacement disorder can be realized in two dif- 
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Fig. 3. Schematic drawing of the structure of CU3AU above and below 
the ordering phase transition at T C =394C°. 



ferent ways. In the first method a randomly selected atom is 
displaced by a random Gaussian distributed amount ('shift', 
Fig.^b). Alternatively the displacement variables associated 
with two different randomly selected atoms are interchanged 
('switch-displacements', Fig. ^Jc). The latter method has the 
advantage that the overall mean-square displacement aver- 
ages for each atom site can be introduced into the starting 
model and these will remain constant throughout the simu- 
lation. 



Occupational disorder: Cu 3 Au 

In this section we will extract chemical ordering information 
from a simulated PDF, which we will refer to as 'data' . This 
allows one to test the capability of the PDF method to extract 
such information for crystalline samples. The basis for these 
simulations is the system CU3AU. The chemical ordering in 
these systems has been studies for quite some time and fea- 
tures as text book example for chemical ordering e.g. in the 
highly recommended book by jCowley (1995 ). The structure 
of Cu 3 Au is cubic, space group Fm3m with a lattice pa- 
rameter of a — 3. 734 A at room temperature. The system 
undergoes a order-disorder phase transition at T C =394C°, 
schematically shown in Figure [|. Above T c the two atom 
types are randomly distributed, below T c copper atoms oc- 
cupy the corners of the unit cell and gold occupies the face 
centers. 



Simulations 



Chemical ordering ordering is most conveniently described 
using the correlation coefficient which is defined as: 



9(1-9) 



(7) 



Pij is the joint probability that both sites i and j are occupied 
by the same atom type and 8 is its overall occupancy. Nega- 
tive values of cy correspond to situations where the two sites 
i and j tend to be occupied by different atom types while 
positive values indicate that sites i and j tend to be occupied 
by the same atom type. A correlation value of zero describes 
a random distribution. The maximum negative value of cy 
for a given concentration 8 is —6/(1 — 8) (P^ = 0), the 
maximum positive value is +1 (Py = 8). Using MC simu- 
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Fig. 4. Change of the calculated PDFs for O13 Au as a function of tem- 
perature. The PDFs were generated from the correlation values listed 
in Table |l| using MC simulations. The difference is plotted below the 
data sets in each panel. Note that the differences for panels (a) to (c) are 
enlarged by a factor of 15, whereas the difference shown in panel (d) is 
not enlarged. 



lations as described earlier, a total of five disordered struc- 
ture were created corresponding to temperatures of: below 
T c (ordered), T=396C°, 420C°, 685C and random distribu- 
tion of Cu and Au. The simulated crystal structures extended 
over 10x10x10 unit cells, containing a total of 4000 atoms. 
The concentration of gold is as expected 8 = 0.25. The tar- 
get correlation values were taken from Chen et al. (1970 ) 
and all eleven neighbor values were used in the simulation. 
The MC simulations were carried out for 200 cycles, i.e. ev- 
ery atom was visited on average 200 times. The target values 
and achieved values after the MC simulations are shown in 
Table |]in columns 'input' and 'MC, respectively. Note that 
the ordered structure of CU3AU was created directly from 
the structural information. The PDFs calculated from these 
structures are shown in Figure f|. The calculation was done 
for x-ray scattering. In order to observe the change in the 
PDF more clearly, PDFs of two subsequent temperatures 
and the resulting difference are plotted as separate panels in 
Figure |j. At high temperature the Cu and Au atoms are ran- 
domly distributed, which corresponds to correlation values 
Cij = 0.0 for all neighbors. At T=685C° we are still well 
above the phase transition at T C =394C° and inspection of 
Table [j] shows a significant negative correlation for the near- 
est neighbor and a positive correlation for the next-nearest 
neighbor. At this temperature the system starts having pre- 
ferred Cu-Au nearest neighbors, however the ordering does 
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Table 1. C113A11 refinement results. 
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not yet extend to larger neighbor distances. The top panel 
in Figure [| shows the PDFs for the random arrangement 
and the situation at T = 685C°. The main differences that 
can be observed for the first few neighbors showing correla- 
tions significantly different from zero. As the temperature is 
lowered towards the phase transition, the differences of the 
PDFs between temperature points become smaller, however, 
as correlations of more distant neighbors deviate from zero, 
changes in the PDF at higher values of r can be seen. Note 
that the difference curves of the top three panels in Figure Q 
are enlarged by a factor of 15. The bottom panel of Figure 
^ show the difference between the PDF just above the phase 
transition and the ordered low temperature state. Here differ- 
ences are clearly visible and this time the difference curve is 
not enlarged. It can be clearly observed that the main change 
of the PDF is not at higher values of r, since in the ordered 
state, non-zero correlations extend over the complete crys- 
tal. 

The aim of these simulation was to explore the capability 
of RMC simulations of these calculated PDFs to recover the 
correct chemical correlation values at least for this idealized 
situation not considering experimental factors. The RMC 
simulation method was described previously. The parame- 
ter a given in equation || was decreased linearly during the 
RMC refinement. This had the effect that initially the ratio 
of accepted 'good' to 'bad' moves was about 2:1 which then 
decreased until no 'bad' moves were accepted in the last cy- 
cles. The refinements were carried out in 'swchem' mode 
which kept the relative abundance of Cu and Au constant at 
its initial values. The refinement was stopped after 15 RMC 
cycles when practically no RMC moves were accepted. 

Results 

The results of the PDF refinements for the ordered struc- 
ture and the three temperature points above T c are shown in 
Figure BL The corresponding correlation values are given in 
Table [yin the column marked 'RMC. Inspection of Figure 
H reveals a nearly prefect agreement between the simulated 
'data' and the RMC result for all data sets except for some 
small differences in case of the ordered low-temperature 
structure. It should be noted, that the differences shown in 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 
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Fig. 5. Result of RMC refinements of Q13 Au. The differences are mag- 
nified by a factor of 5 and plotted below the data in each panel. For 
details see text. 



Figure ^ are enlarged by a factor of 5. When comparing the 
resulting correlation values, one finds a very good agreement 
between the initial correlations (column 'MC') and the re- 
finement results (column 'RMC'). The most noticeable dif- 
ferences again are for the ordered structure. It is, however, 
not surprising that the RMC refinement is not able to pro- 
duce the ordered structure exactly without some residual de- 
viations. The results clearly show that chemical ordering in- 
formation from crystalline materials can be obtained using 
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2 3 4 5 6 7 8 9 10 11 12 13 14 
r(A) 



Fig. 6. PDFs of simulated structures: (a) occupational disorder only 
and (b) occupational and displacive disorder. In both panels the PDF 
of the undistorted structure is shown for reference. The differences are 
plotted below the data in each panel. For details see text. 
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Fig. 7. Results of RMC refinements: (a) occupational disorder only and 
(b) occupational and displacive disorder. The differences between are 
magnified by a factor of 5 and are plotted below the data in each panel. 
For details see text. 



RMC refinements of PDF data. However, one should keep in 
mind, that these simulations were carried out to investigate 
the principle possibilities of this technique, excluding prob- 
lems introduced by noise and systematic errors one would 
encounter when using real experimental data. 



Occupational and displacive disorder 



Table 2. Results of RMC refinements of disordered structures showing 
occupational as well as displacive disorder. The bond lengths are given 
in A. 
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The next step is to add distortions to a defect structure show- 
ing chemical short range order. In this example, the simu- 
lated test structure had no 'real' analog. For convenience, 
we used the same concentration of Cu and Au atoms as in 
Cu 3 Au. 



Simulations 

The disordered structures were created in two steps: First 
the desired chemical ordering was introduced the same way 
as described in the previous section. As a second step, the 
atoms were displaced using MC simulations and a Hooke's 
law potential as described above. The starting structure used 
is cubic primitive with a lattice parameter of a = 3.75A. The 
model crystal extends over 32x32x32 unit cells giving a to- 
tal of 32768 atoms. The site at (0,0,0) is randomly occupied 
with 75% copper atoms and 25% gold atoms. The chemical 
arrangement is characterized by a negative nearest neighbor 
correlation of cioo = —0.3 and a positive second nearest 
neighbor correlation of cuo = 0.5. In other words nearest 
neighbor Cu-Cu and Au-Au pairs are avoided whereas sec- 
ond nearest neighbor pairs of the same atom type are favored 
compared to a random arrangement. The correlation coeffi- 
cients actually achieved by the MC simulation are listed in 
Table |[ The PDF calculated from the structure containing 
only the occupational short-range order is shown in Figure 
^a) superimposed on the PDF of the random structure. The 
difference between both data sets is plotted as well. 

The effect of thermal motion on the PDF is simulated 
by convoluting with a Gaussian. In addition the atoms are 



shifted from their ideal posi tion by the MC sim ulation in- 
troducing a size-effect type (Butler et al., 1992) distortion. 



In this case the resulting nearest neighbor distances are is 
dcu-au = 3.73A and dcu-cu,au-au = 3.77A (Table ||). Note 
that the average distance is the lattice parameter a = 3. 75 A. 
A total of 200 MC cycles were computed to create this struc- 
ture. The calculated PDF is shown in Figure ||in the bottom 
panel. The PDF of the undistorted starting structure is plot- 
ted in the same panel and the difference between both is 
shown below the data. 

Two different RMC refinements were carried out: Run O 
using the 'data' containing chemical ordering only and Run 
D using the 'data' containing additional displacements. Run 
O was carried out as a reference of the chemical correla- 
tions that can be refined when no displacements are present. 
The refinement was carried out exactly as described in the 
last section. The resulting PDFs are shown in Figure ^ and 
the achieved correlation values for run O are listed in Ta- 
ble^. As for CU3AU, we observe an excellent agreement be- 
tween 'data' and the calculation. The interesting question is 
now, can one extract the chemical short range order equally 
well with the size-effect distortion present, and can one ex- 
tract the magnitude of those distortions. In order to model 
occupational and displacive disorder simultaneously during 
the RMC refinement, we employed the following strategy 
which was already found to b e successful when refining sin - 
gle crystal diffuse scattering ( feoff en and Welberry, 1997b ). 
The occupational shifts and displacement shifts were carried 
out alternately with no more than 10% of the crystal sites 
being visited before switching between the two modes. As 
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before, the parameter a was decreased during the RMC re- 
finement in such a way that initially a significant amount of 
'bad' moves was accepted, however, in the final few cycles, 
no moves worsening the agreement between observed and 
calculated PDFs are accepted. The simulations were carried 
out for a total of 15 cycles when the agreement did not sig- 
nificantly improve any further. The resulting PDFs of run D 
are shown in Figure ^ and the achieved correlation values 
are listed in Table 0. 



Results 

Inspection of Figure ^a shows a very good agreement be- 
tween the RMC result and the 'data' for run O. The chemi- 
cal correlations of the resulting structure match the expected 
ones nearly perfectly (see Table ^) similar to the result found 
for the CU3AU simulations presented in the last section. The 
focus in this part is on run D where chemical and displacive 
contributions were present. The refined PDF shows a good 
agreement with the 'data' (Fig. ^|b), the only noticeable dif- 
ferences are at the nearest neighbor distance and at twice 
that distance. This indicates some residual differences in the 
exact distribution of distortions within the crystal. It is in- 
teresting to note that the successful refinement of a system 
showing displacements required a much larger model con- 
taining in excess of 32000 atoms compared to only 4000 
used for the simulations of CU3AU. In addition to the static 
displacements introduced as size-effect, the contribution of 
thermal vibrations to the PDF needs to be accounted for. 
Various test runs showed that modeling the thermal com- 
ponent via convolution of the PDF with a Gaussian gives 
better results compared to randomly displacing all atoms in 
the model structure. Inspection of the achieved correlations 
and displacements listed in Table || shows an overall good 
agreement with the expected values, although the nearest 
neighbor correlation does not match as good as for run D 
compared to the result of run O. The resulting distortions 
match the expected values as well. In summary, the chemical 
correlations as well as magnitudes of distortions were suc- 
cessfully extracted from the PDF using RMC simulations. 



Conclusions and outlook 

Two important conclusions can be made from the results 
presented in this paper: First the chemical short range or- 
der information can be extracted from the PDF using RMC 
simulations even in cases where additional displacive disor- 
der is present. Furthermore, one can also extract the mag- 
nitude of the specific distortions very reliably. For this sim- 
ple isotropic disorder, the extracted correlations for the un- 
derlying chemical ordering from the PDF are even slightly 
closer to the expected values than in a similar set of RMC 
refinements of single crystal diffuse scattering (Proffen and 
Welberry, 1997b). It is interesting to compare the diffuse 
scattering calculated from the input and resulting structure 
of run D. The calculated diffuse scattering for the zeroth 
layer of the hk-plane in reciprocal space is shown in Fig- 
ure |[ The patterns show a qualitative agreement between 
'data' and refinement result, which strengthens the results 





12 3 4 
[h 0] 



Fig. 8. Diffuse scattering from simulated structures: (a) input structure 
for simulation 'Both', (b) corresponding RMC result. Units are recipro- 
cal lattice parameters. The calculation was done using neutron scatter- 
ing. For details see text. 



obtained from refining the PDFs. In case of refining the dif- 
fuse scattering directly one obtains a much better mat ch with 
the observed data (see Proffen and Welberry, 1997b). How- 
ever, although the diffuse scattering patterns shown in Fig- 
ure H are in general similar, closer inspection reveals an in- 
teresting difference: The pattern corresponding to the input 
structure has in general sharper more localized diffuse scat- 
tering whereas the pattern calculated from the RMC result is 
somewhat broader. This indicated that the defects in the re- 
sulting structure are not correlated over long enough lengths. 
Although the PDF contains diffuse scattering information as 
well, there is a certain complementarity between single crys- 
tal diffuse scattering and PDF data. We are currently inves- 
tigat ing this relationship on the system Ini_ z Ga-,;As (J eong 



et al.). This opens the possibility of a joined refinement of 
single crystal diffuse scattering and PDF data obtained from 
powder data. 

The 'data' used for the refinements in this paper were cal- 
culated from given disordered structures, which allowed us 
to judge the results since the 'correct' result was known. Ad- 
ditionally those data were free from experimental errors and 
noise, which of course will make refining of PDFs obtained 
from real experimental data more challenging. The simula- 
tions presented here were carried out using the total PDF, 
i.e. all atom pairs are included in the PDF. Using anomalous 
scattering or by combining neutron and x-ray data, one can 
separate out the different partial PDFs. For a simple binary 
compound consisting of atoms A and B, the corresponding 
partials would contain only AA, BB or AB (same as BA) 
atoms. More details can be found in a recent study of the 
total and differential PDF of Ini-^Ga^As by Petkov et al. 
(1999). The program DISCUS used here can also be used to 
calculate or refine any partial PDF of a given structure. A 
different definition for partial PDFs was proposed by B hatia 
and Thornton (1970). They separated the different terms in 
a binary case into a topological term (the usual PDF), the 
radial concentration correlation function (RCF) containing 
information about chemical fluctuations and a cross-term. 
Current efforts are to include the calculation and refinement 
of the RCF into the DISCUS program. 

In summary, we have shown in this paper that RMC re- 
finements of the PDF of disordered crystalline materials can 
extract chemical short-range order information for systems 
showing only occupational disorder as well as for systems 
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showing occupational and displacive disorder in combina- 
tion. The usage of faster and usually simpler powder diffrac- 
tion measurements allows one to extract this information 
more quickly and cover e.g. many temperature points around 
a phase transition. However, for certain problems, one will 
still need to use single crystal diffuse scattering to deter- 
mine the details of the defect structure or even perform a 
combined refinement of PDF and single crystal data. Mea- 
surements of the system CU3AU using synchrotron radiation 
are planned and will be the first 'real' test using the PDF 
method to extract chemical short range order parameters for 
this crystalline system. 
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